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ABSTRACT 

We measure the decay rate of the mid-IR himinosity from type la supernova 201 Ife between six 
' months and one year after explosion using Spitzer/IRAC observations. The fading in the 3.6 /im 

\ channel is 1.48±0.02 mag/100 , which is similar to that seen in blue optical bands. The supernova 

, brightness fades at 0.78±0.02 mag/100'^ in the 4.5 /im channel which is close to that observed in the 

' near-IR. We argue that the difference is a result of doubly ionized iron-peak elements dominating the 

I bluer IRAC band while singly ionized species are controlling the longer wavelength channel. To test 

. this, we use Large Binocular Telescope spectra taken during the same phases to show that doubly 

[•-[- ' ionized emission lines do fade more slowly than their singly ionized cousins. We also find that [Co III] 

' emission fades at more than twice the radioactive decay rate due to the combination of decreasing 

excitation in the nebula, recombination and cobalt decaying to iron. The nebular emission velocities 
' of [Fe III] and [Co III] lines show a smaller blue-shift than emission from singly ionized atoms. The 

Si II velocity gradient near maximum light combined with our nebular velocity measurements suggest 
SN 201 Ifc was a typical member of the 'low velocity gradient' class of type la. Analyzing IRAC 
photometry from other supernovae we find that mid-IR color of type la events is correlated with the 
' early light curve width and can be used as an indicator of the radioactive nickel yield. 

Subject headings: supernovae: general — supernovae: individual (SN 2011fe, SN 2009ig, SN 2008Q, 
SN 2005df 
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1. INTRODUCTION 



C§ ' Typ e-la supernovae (SNIa) are widely thought to be the thermonuclear explosions of degenerate white dwarfs (see 
lLivioll200Q f or a r eview). Their hi gh luminosity and uniformity make SNIa good cosmological distance indicators 
( Riess et all 119981 : iPerlmutter et al.] |l999il , but a conclusive theoretical model of their explosion mechanism remains 



elusive. For instance, it is an open question on whether the progenitor consists of a single degenerate white dwarf or 
two (or both). 

'"J We do know that SNIa are bright because they synthesize a sizable amount of radioactive ^^Ni. Typically, SNIa 
produce half a Solar mass of radioactive nickel which decays to cobalt (^1 week time sc ale) that then decays to stable 
^^Fe (~three month time scale) emitting energetic 7-rays and positrons ()Arnettl [l982|) . Early spectroscopy shows a 
significant amount of intermediate mass eleme nts (IME) in the outer layers while late-time spectra are dominated by 
CN singly and doubly ionized iron-peak elements (jFilippenkol 1 1 99 7|) . 

' SNIa explosions may begin with a subsonic burning front, creating a deflagration of the carbon/oxygen white dwarf, 
■ but a deflagration alone yields too little energy. Alternatively, a pure detonation would create too much radioactive 
, yield and less IME than t ypically observed. A "delayed detonation" (DD) better matches the observations. In the DD 
' scenario (iKhokhlovlll99l . fusion be gins as a deflagration but switches to supersonic burning after some of the white 
dwarf has expanded, allowing carbon fusion at low density. 
Recent work has suggested th at the de nse defiagration region can be offset in position and velocity from the larger, 
?-j ■ lower density detonation region (Maeda et al. 2010a). This proposition explains a number of observed SNIa properties 
near maximum and at later phases and that some diversity in SNIa may simply be due to random viewing angles. 

At late times, when the synthesized material becomes optically thin, the light from SNIa continues to be powered by 
radioactive decay. Months to years after the explosion, the long-lived ^^Co isotope keeps the nebula glowing by energy 
deposited from positrons while the 7-rays escape without depositing energy. The photometric decay is expec ted to be 
exponential with a e-folding time of 111.3 days or 0.98 mag/100'^, matching that of ^^Co deca v. ButlLair et al. (2006) 
found the optical decline in the B and V bands is significantly faster than ^^Co . In contr ast, SoU erman et al.i (i200j) 
showed the near-IR bands fade very slowly a year after explosion. However, Solle rman et a l. (2004) also demonstrated 
that the bolometric luminosity, or the total fiux across all bands from the optical to the IR, does decay at the expected 

^^Co energy deposition rate . 

SN 2011fe ()Nugent et all 1201 ID was the closest SNIa in 25 years (/i « 29.04; iShappee fc Stanekl [20T1I) and was 
extremely well studied at early epochs. It appears to have been a normal SNIa in terms of light curve, luminosity and 
spectral development. Here, we present Spitzer IRAC photometry of SN 2011fe between six months and one year after 
optical maximum. The brightness of SN 201 Ife created a unique opportunity for detailed study of a type la event 
in the mid-IR, but it also appeared close to the Earth's ecliptic pole, allowing several months of continuous Spitzer 
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visibility. To put these results into context, we also compare Spitzer observations of SN 201 Ife with data for other 
SNIa. 

2. OBSERVATIONS 
2.1. Spitzer /IR AC 

As part of a program to study the late-time properties of SNIa, we requested low-impact target-of-opport unity obser- 
vatio ns of SN 2011fe with the Spitzer Space Telescope using the Infrared Array Camera (hereafter IRAC: IFazio et al] 
I2004D . IRAC photometry was also obtained through Spitzer program 80196. The location of the supernova was first 
visible to Spitzer at the end of January 2012 and remained visible until the following September. Four visits were exe- 
cuted over the visibility period providing the best light curve of any SNIa in the mid-IR. The supernova was observed 
in IRAC Channel 1 (CHI) and Channel 2 (CII2) which have effective wavelengths of 3.6/im and 4.5/.tm respectively. 

Our program also obtained IRAC photo metry of SN 2008Q a nd SN 2009ig, and we re-reduced archival IRAC 
data on SN 2005df previously pres ented bv iGerardv et all (I2007D. These three SNIa were well-observe d at optical 
wavelengths around maximum light (|McClelland et al. Ilin preparation! : iFo lev et al.ll20l3lMilne et al.ll2010l ) and provide 
good comparisons to SN 201 Ife. 

We downloaded post-BCD images through the Spitzer Heritage ArchiveQ. Each SNIa was identified in the IRAC 
images by its celestial coordinates^- Aperture photometry was performed after the pixel values were converted to /iJy. 
For SN 201 Ife, we used the recommend 5-pixel aperture radius with annular sky radii of 12 to 20 pixels. Pre-explosion 
Spitzer observations of host galaxy MlOl conducted in 2004 provided background estimates for both channels. The 
position of SN 2009ig in host galaxy NGC 1015 warranted a smaller 3-pixel aperture radius, while sky backgrounds 
were calculated via averaging the brightness fluctuations of multiple random apertures over areas without obvious 
point sources. The same was done for SN 2008Q in host galaxy NGC 524 and SN 2005df in NGC 1559. Aperture 
corrections were taken from the IRAC Instrument Handboolfl and the results are presented in Table [T] 

2.2. Large Binocular Telescope 

Spectr a of SN 20 life wer e obtained with the Large Binocular Tel escope (LBT) and M ulti- Object Dual Spectrograph 
(MODS; i Pogge et al.ll2010l ) on several epochs during its decline (jShappee et al.l[201^ . Three of the spectra, 2012 
March 24, April 27 and June 12 (UT) were obtained over the span of the Spitzer observations and we use these to 
estimate the decay rates of SN 2011fe at optical/near IR wavelengths. The epochs correspond to ages after S-band 
peak light of 196, 230 and 276 days respectively. 

The spectra were obtained through a 1.0 arcsecond wide slit that was oriented to the parallactic angle for the March 
observation, taken at an airmass of 1.15. The April and June observations were obtained at airmasses less than 1.1 so 
that differential slit losses were minimized. 

The spectra were wavelength-calibrated using Neon and Argon emission line lamps taken during the June run. The 
wavelength solution was then shifted based on the sky line centroids. The MODS spectrograph is very stable and the 
wavelength shifts applied were less than 2 A. 

MODS flux sensitivity is also very stable, so no flux calibration was applied. Instead, we directly compared ex- 
tracted electron count rates (proportional to the 'analog-to-digital units') to avoid further uncertainties caused by flux 
calibrations using different standard stars under varying conditions. To test the stability of the sensitivity, we divided 
MODS spectra of the standard star HZ 44 obtained in April and June and found the fluxes differed by less than 5% 
over the spectral range between 3500 A and 8500 A. We therefore expect systematic uncertainty on the pixel by pixel 
decay rates to by less than 0.05 mag/lOO"^. At wavelengths shorter than 3500 A and longer than 8500 A the difference 
between the standard star spectra reached 10%, so the decay rates are not as reliable at the extreme ends of the 
spectra. 

The MODS spectra of SN 2011fe are shown in Figure [TJ In general, the spectra are typical of a type la supernova a 
couple hundred days after explosion. Between 4000 A and 5500 A there is a clump of emission dominated by forbidden 
lines of singly and doubly ionized iron. Weaker lines from cobalt, iron and nickel are clearly seen long ward of 5500 A. 
These are some of the highest signal-to- noise (S/N) spectra of a type la supernova at late p hases ever obtained and 
they already been used to search for hydrogen emission from SN 2011fe (|Shappee et al.l[20m . 

The relative count rates across the spectra are stable from epoch to epoch, but the total efficiency depends on seeing 
and cloud cover and these were not well controlled. To simulate a broadband filter, we multiplied the spectra with a 
Bessell i3-band filter function and integrated the result. From these values we adjuste d the relativ e ADU to match 
the S-band decay rate of 1.32 ± 0.02 mag/lOO'^ found from the late-time photometry of iMunari et al. (2013). 

LBT spectra of SN 2011fe were also obtained on 2012 May 1 (UT, age of 234 days) with the LUCIFER infrared 
spectrograph. The spectra covered the J and H bands using three grating settings. Only one epoch of near-IR spectra 
was obtained so the decay rates at these wavelengths could not be estimated, but the velocity offsets of the two strong 
[Fe II] lines are measured. 

3. ANALYSIS 

The mid-IR photometry of SN 201 Ife starting at an age of five months and ending about one year after explosion 
is shown in Figure [2j The light curves are consistent with an exponential decay, so we fit them with a function of 

http:/sha. ipac.caltech.edu/applications/Spitzer/SHA 
^ littp://ned. ipac.caltech.edu/ 

® http:/ /irsa.ipac. caltech.edu/data/SPITZER/docs/irac/iracinstrumenthandbook/ 
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Figure 1. LBT/MODS spectra of SN 2011fe obtained during tlie 2012 Spitzer observations. The spectra are not flux calibrated, but are 
plotted as the logarithm of the raw counts per pixel. The MODS dichroic divides the red and blue sides of the spectrograph at 5700 A and 
this is shown as a gap in the plot. Some emission features are seen to fade at a different rate than others. For example, the emission peaks 
around 6000 A decline significantly faster than the features at 7200 A. 

the form J- = ^2306"'*"^^*^-'/'^, where T is the observed time varying flux, t is the age since B maximum, A230 is the 
flux at the fiducial time of 230 days, and r is the e-folding time of the decay. The decay rate can also be described in 
magnitudes per 100 days by A = 108.57/t. 

The resulting measured decay rates for SN 2011fe are given in Table [2j The 1.48±0.02 mag/lOO'^ seen in CHI is 
significantly steeper than the 0.78±0.2 mag/lOO'^ measured for CH2. The decay observed in the two Spitzer bands 
straddle the '^^Co rate of 0.98 mag/lOO*^. The decay rate for CHI is similar to that typically seen in B and V bands, 
whil e CH2 has a decay e ven more shallow than the / band (0.99 ± 0.04 mag/lOO'^; via least-squares fit to photometry 
from iMunari et all I2013D . 

The long cadence and greater distance to the other Spitzer-observed type la events resulted in larger uncertainty 
in the measured decay rates compared with the SN 2011fc data. The CHI decay rates for SN 2008Q, SN 2009ig and 
SN 2005df all appear slightly slower. All the supernovae are consistent in that the CHI rate is steeper than the CH2 
decline rate. 

We created a color index based on the fiux ratio between the two IRAC channels, or 

CHI - CH2 = -2.5 log f ) + {ZPcYii - ZPcm) mag, 

\-rcm J 

where -Fchi,ch2 are the fluxes measured in /iJy in CHI and CH2, and ^Pchi,ch2 are the magnitude zero points for 
those passba nds, listed in the IRA C Instrument Handbook. Extinction differences between these filters is assumed to 
be negligible ([Flahertv et al.|[2007D . Because the IRAC decay rates are so different, the CH1-CH2 color varies rapidly 
in time so we have chosen a fiducial time of 230 days to estimate the color parameter. This age was chosen because 
the Spitzer archive contains several SNIa observed near this age, thus reducing any errors caused by interpolation. 

The estimated color at 230 days shows a wide range among the supernovae in our sample. To better understand 
the cause of this diversity, we sought to compare the 230 day color wi th early-time ligh t curve shape, namely the 
parameter Ami5(i?j3. For SN 2009ig we use the ATOi5(i?) measured bv lFolev et al.l (|2012| ) and for SN 2011fe we use 



^ Ami5(_B) HPhillipslll9931 ) is defined as the change in B magnitude from maximum to 15 days post maximum. It is related to the peak 
luminosity of a SN la through the Phillips Relation. 
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Figure 2. 2011fe IR light curve computed from IRAC ptiotometry. Tlie Ctiannel 1 decay slope is consistent with typical SN B and V 
decay slopes, while Channel 2 decays faster than the ^'^Co energy deposition rate. If the decay rate in each filter remains constant, 2011fe 
should see equal brightness in both channels at ~ 1 year after maximum. 



the value from iVinko et all (120121). Both SN 200Qdf and 2008Q have excellent early light curves (|Milne et al.l 120101 : 
IMcClelland et al. Ilin pr eparation^ and we measure their Ami5(_B) directly. The Ami5(_B) values and the measured 
IRAC colors are given in Table [2j 

The four supernovae, 2005df, 2008Q, 2009ig and 2011fe, show a wide range of early light curve decline rates, and so, a 
wide range of synthesized radioactive nickel masses. We plot our measured IRAC colors at 230 days versus Ar7ii5(_B) in 
Figure |3l There appears to be a robust correlation (albeit with only four data points) indicating that the CHI to CII2 
flux ratio decreases with increasing Arni5(_B) parameter. Alternatively, we see that the CHI flux increases relative to 
CH2 as the amount of synthesized radioactive nickel increases since Ami5{B) correlates with nickel yield. Based on 
this handful of supernovae, the IRAC color appears to be a good indicator of peak optical luminosity, although this is 
unlikely to be useful for cosmological measurements. 
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Figure 3. Relation between late-time IRAC color and early-time decline rate. Magnitudes are calculated according to the photometric 
zero points given in the IRAC Instrument Handbook. The dashed line represents the least-squares linear fit to this relation. 
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4. DISCUSSION 
4.1. IRAC Decay Rates 

The radically different late-time decay rates observed at 3.6/im an d 4.5^tm is puzzling given that total flux fro m the 
optical to near-IR appears to match the ^^Co energy injection rate (jSollerman et al.ll2004l : iLeloudas et al.ll2Q09f ). The 
effect is also seen when comparing the faster decay in the B and V bands to the modest /-band decay rate between 
six months and a year after explosion. The fast optical decline is likely compensated by a slow near-IR decay which 
after a year becomes a spectacular rise in the importance of the J and H band fluxes. 

What do the fast declining wavelength bands {B, V, CHI) have in common and how are they different from the 
slower (/, J, H, CH2) fading bands? As most of the SNIa spectrum at late-times, the IRAC bands are expected 
to be dominated by emission from iron-peak elements with a mixture of singly-ionized and doubly ionized species. 
iGerardv et al.l (2007) predicted that spectrum covered by CHI would show strong [Fe III] 3.229 /zm and [Co III] 3.492 
lim. In contrast, Macda et al. (2010a) predicts several strong [Fe II] should dominate the wavelength region spanned 
by CH2. 

We propose that the difference in the observed IRAC decay rates is a result of doubly ionized species being the 
dominate light source in CHI and singly ionized iron-peak species supplying the majority of flux in CH2. The relative 
amounts of emission from singly and doubly ionized species should also vary in the shorter wavelength bands, and 
we note that the strong [Fe III] 4700A blend straddles the B and V band transition and contributes flux to both 
fast-decaying bands. In contrast, the slow-fading /, J and H bands are predominately fed by singly ionized atoms. 
We predict the if-band around 2.2/im should fade like IRAC CHI since it contains a strong [Fe III] line at 2.218/im. 

Si ngly and doubly ionized i ron-peak elements exist in very different regions of SNla nebulae. Using the W7 prescrip- 
tion (jThielemann et a"Llll986D . one-dimensional SNIa nebular models agree that [Fe II] (and [Co II] which has a similar 
ionization energy) atoms dominate i n the inner dense regions while doubly ionized iron occupies the outer volume (e.g. 
iMazzali et al.l[20lTI : iLiu et aLlll997| ). Generally, the dense regions cool quickly and doubly ionized species reco mbine 
while cooling in low density regions is inefficient and they remain dominated by doubly ionized gas. SoUcrman et al.l 
(|2004f l however, notes that photoionization has a major effect on the ionization structure and it is not well modeled 
i n the nebular spe c tra (|Kozma fc Franssonlfl99l . 

iSollerman et al.l (|2004| ) modeled SNIa nebular spectra with and without accounting for photoionization. Without 
including photoionization in the calculation, the [Fe III]/[Fe II] ratio was too small compared with the real spectrum of 
SN 2000cx at one year old. When including photoionization, the model predicted too much [Fe III] which suggests that 
real supernovae are intermediate between these two extremes. It is possible that the effects of photoionization in SNIa 
nebulae decrease with time as the expansion reduces the density. This should reduce the importance of the doubly 
ionized species through recombination while increasing the flux of singly ionized species. As the nebula transitions 
away from significant photoionization, we expect the filters dominated by lines from doubly ionized species will fade 
more quickly than ^^Co. Meanwhile, the bands dominated by singly ionized atoms will fade more slowly than A = 0.98 
mag/lOO*^ with the difference in decay rates a clue to the recombination rate. 

If the CHI decay rate is a result of both the excitation decreasing at the ^^Co rate and the recombination rate from 
[Fe III] [Fe II] , then we can approximate it by 

exp{-t/Tcm) = exp{-t/Tco) ■ exp{-t / t^^c) , 

where rcHi7 tco and t^cc are respectively the e- folding times of CHI, ^^Co and [Fe III] [Fe II] recombination. Assum- 
ing coUisional ionization are negligible at late times and photoionization is becoming unimportant, the recombination 
should go approximately as exp^—ncant), where rig is the electron number density and an is the recombination rate 
coefficient. The recombinati on coeflicient is an « 3 x 10""'^^ cm'^ s~^ and is nearly constant for temperatures between 
10^ and 10^ °K (|NahaHll997l ). We can solve for the electron density, = l/aj?[l/TcHi - 1/tco] and find = 1.8 x 10^ 
cm~^. Assuming most of the iron peak atoms are doubly ionized, the mass density corresponding to this recombination 
rate is then approximately 10^^^ g cm^"^. 

We can compare this recom bination density to the density of the standard W7 model scaled to 250 days after 
explosion (jMazzali et al.lf2"011l ). The scaled W7 density ranges between 10~^^ g cm~'^ in the inner regions to 10^^*" 
g cm~^ out at 10000 km s~^. Considering the assumptions made in the calculation, this is surprisingly close to the 
density for a pure recombination process removing doubly ionized species from the nebula. 

4.2. Optical Decay Rates from LBT Spectra 

The ratio of the June 12 and March 24 count rates as a function of wavelength is shown in Figure |4] in units of 
mag/lOO'^. The ratio has been scaled to match the B-band decay rate discussed in Section [51 Clearly, the decay rates 
vary by spectral region and by specific feature. Some of the small-scale variation is due to shifts in emission centroids 
with time. For example, the bright blend of [Fe III] lines at 4700 A shifts by about 10 A redward over the 80 days 
between spectra. This is shown as a decay ra te of 1.4 mag/ 100'^ o n the blue side of the feature and a 1.1 mag/100'^ 
decay on the red side. This drift was noted bv'M aeda et al.l (|2010a| ) and appears typical in SNIa. The total fiux of the 
feature fades by 1.41 mag/lOO'^ which is not surprising given that we have scaled the S-band fiux to 1.32 mag/100'^ 
and the 4700 A line contributed significantly to this filter. 

In contrast, the weak features around 6000 A faded at twice the ^^Co decay rate and show no shift in the line 
centroid. The 5900 and 6160 A lines have been attributed to [Co III] emission, Wt the 5900 A line may also contain 
some Na 1 D fiux. 
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Figure 4. Decline rate A in units of mag/100'^ from the LBT/MODS spectra from 190 days to 270 days past maximum. Tfie optical 
spectrum of 2011fe is included for reference. The horizontal dotted line at 1.41 mag/IOC* represents the calculated _B-band decline rate, 
while the horizontal line at 0.98 mag/IOC^ is the energy deposition rate, equal to ■"'"Co decay. We note that the large fading rate seen from 
5700 A to 6400 A is Co III which has an additional loss channel through radioactive decay to iron. 

As noted in Figure [U the blend of bright lines around 7200 A fade rather slowly compared with the [Fe III] and 
[Co III] lines in the blue. The emission around 7200 A is often attributed to [Fe II] and [Ni II] lines. The line at 8600 A 
also shows a slow decay rate and it too has been identified as an [Fe II] feature. 

In general, the singly ionized species show the slow decay rates while emission from doubly ionized atoms decays 
significantly faster than the radioactive decay rate (see Figure [5]) . This supports our interpretation of the decay rates 
we observed in the Spitzer bands and that recombination is the likely sink for the doubly ionized species. 

[Co III] emission is hit with a triple whammy. The overall energy available for line emission decreases at the ^^Co 
radioactive decay rate. In addition, recombination of doubly ionized atoms reduces their numbers as demonstrated 
in the decay rate of [Fe III] lines. But much of the cobalt is radioactive, so its numbers are decreased at twice the 
radioactive decay rate plus the recombination loss rate. 

The emission at 6550 A has been identified as from [Co III], but it appears to fade at a slower rate the the other 
[Co III] lines. The very strong shift in peak wavelength as the line fades suggests it is a blend with a singly ionized 
emission feature and that combination makes it appear to fade more slowly than the other [Co III] lines. T he NIST 
Atomic Spectra Databas43 points to strong [Fe II] emission at 6446, 6456 and 6516 A, while iLiu etaH ([1991 suggests 
a combination of [Fe II] and [Ni II] in addition to [Co III]. This feature's profile is flat-topped in the June spectrum 
and this makes it difficult to obtain a precise peak velocity. 

4.3. IRAC Color Variation in Type la Supernovae 

Our supposition that emission from doubly ionized species decays at a faster rate than those from singly ionized 
atoms appears to be supported by the optical spectr a. Since IRAC CHI s hows a rapid decay rate, then doubly ionized 
species should be dominant around 3.6 /im. Indeed. fCerardv et al.l ()2007[ ) suggested Fe III and Co III lines contribute 
significantly to the CHI band. Thus, the CHI— CH2 color observed in the SNIa might be a good indicator of the 
ionization ratio within a SNIa nebula. 

We can ask if the correlation between the IRAC color and Ami5{B) seen in Figure p implies that a [Fe III] / [Fe II] 
ratio correlates with nickel yield. SN 2003hv possessed notably strong [Fe III] 4700 A at late times (|Mazzali et al.l 



® http:/ /physics. nist.gov/PhysRefData/ASD/lines_form. html 
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Figure 5. Decay rate of the total flux for individual emission lines, measured between 24 March, 2012 and 12 June, 2012. The approximate 
BVRI (Munari ct al. 2013), CHI and CH2 decay rates of SN 2011fe are given for reference. 



1201 If ), though the SN had a large l\rai^{B) of 1.6 (iLeloudas et al.|[2OO90 . so probably synthesized a relatively small 
amount of radioacti ve nickel. This trend is the opposite of that implied by [Fe III] dominating the CHI bandpass. 
iMazzali et al.l ()201lD argued that the very strong [Fe III] emission seen in SN 20031iv might be due to an unusually 
low nebular density and that the explosion might be peculiar. The [Fe III]/[Fe II] ratio is probably a good indicator of 
nebular density because the recombination rate is so sensitive to density. But it is not clear why the nebular density 
at a fixed age should be related to Ami5(i?). 

If the fast decaying [Co III] 3.492 /im line dominates the CHI bandpass with contributions from slower decaying 
emission such as [Fe II] and [Fe III], then the flux would appear to decay at an intermediate rate over our Spitzer 
observation window. Such a scenario would explain Figure [3] because the CHI— CH2 color would correlate with the 
[Co III]/[Fe II] ratio and that should be a good indicator of the radioactive nickel yield. This supposition awaits L-band 
spectra to identify the emission features in the CHI bandpass. 

We note that the IRAC color could be matched in the optical using the SDSS r — i color index for SNIa at late 
times. The SDSS-r band encompasses fairly uncontaminated Co III features while SDSS-i band mainly sees emission 
by Fe II. Such a color at a fixed age should correlate with peak luminosity. 

4.4. Nebular Velocity 

iMaeda et al.l (|2010aD suggested that asymmetry in SNIa explosions could be ma nifested in the ne b ular spectra. 
While most SNIa appear to show little asymmetry in the outer layers of the ejecta (jWang et al.l 119961 ). IMaeda et al.l 
(l2010b( ) showed that differences in velocity shifts between ions associated with dense, inner layers could be quite large 
when compared to velocity shifts of the outer layers. The theory they develop is that of a dense, inner region rich 
in deflagration products that is offset in position and velocity from the low-density, symmetric outer region. In this 
scenario, doubly ionized iron-group lines tend to trace outer low-density regions, while their singly-ionized variants 
trace the dense inner regions. 

We measured the wavelength centroids of several major features of the LBT/MODS spectra and found two distinct 
groups: [Fe III] and [Co III] lines tended to be blueshifted by ^400 km s~^, while [Fe II] and [Ni II] features were 
blueshifted between 1000 and 1500 km s^^ . A list of these measured lines is given in Tabled Features made from 
blends of singly and doubly ionized lines (e.g. the ~5000 and ^5250 A [Fe II] -I- [Fe III]) appear to show higher 
blueshifts, as well, suggesting that their singly-ionized components are dominant. The highly blended ~6550 feature 
proved difficult to mode l in this man n er; th e lack of a narrow, discernible peak meant a large error in identifying 
individual lines, though iBowe rs et al.l (|1997f ) identified [Co III] 6578 A as an important component, and the NIST 
Atomic Spectra Database points to major [Fe II] contributions here, as well. Synthetic spectra and improved atomic 
transition information is necessary to better determine the shift of individual lines. 

Despite the problems of line blending in the optical, it is clear that the velocity shifts of the [Fe III] 4700 A line 
and the [Co III] features at ~5900 and ^^6200 A are significantly smaller than the offsets of more than —1000 km s~^ 
observed in the singly ionized emission lines. Both [Fe III] and [Co III] emission likely originate in the the outer, 
symmetric regions of the nebula. 
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The large difference in velocity between the [Fe II] and [Fe III] features places SN 2011fe firmly in the "low-velocity 
gradient" group (see Figure |6]). We confirmed this by examining the nebular velocity (defined bv lMaeda et al.|[2010al 
as the average of the [Fe II] 7155 A and the [Ni II] 7 378 A velocity shifts) and the early-time Si II velocity gradient, 
which we calculated using the spectra presented by iParrent et alJ (j2012[ ) and made publicly available through the 
Weizmann Interactive Supernova data REPositorjIj. The observed velocity shift of the s i ngly io nized nebular lines 
implies a viewing angle of roughly 70° from the symmetry axis as defined bv lMaeda et"all (|2010bl ). 




± 



-3000-2000-1000 1000 2000 3000 
v„ei, [km/s] 

Figure 6. Comparison of nebular velocity with respect to host galaxy rest frame to the early-time velocity gradient of Si II (see lMaeda et ahl 
I2010bl for details). SN 2011fe is labeled in blue amongst the green Low Velocity Gradient (LVG) SN. The gray, "Intermediate" SN include 
faint events SN 1986G and SN 2007on, while the High Velocity Gradient (HVG) is colored pink. Peculiar SN such as SN 2004dt have not 
been included in this plot. 



5. CONCLUSIONS 

We have analyzed mid-IR photometry and optical spectroscopy of the nearby type la supernova 201 Ife at ages 
between six months to a year after the explosion. We found the following: 

• The luminosity in the Spitzer/IRAC bands are well fit by an exponential fading. The observed decay rate was 
1.48±0.02 mag/lOO'^ in the 3.6 /xm band and 0.78±0.02 mag/lOO'^ at 4.5 /xm. The fading rate in the CHI fiher 
is similar to that seen for many SNIa in the optical B and V filters. The CH2 decline rate is slower than the 
^^Co radioactive decay rate. 

• Using SN 2011fe and three other SNIa observed by Sptizer, we found that the ratio of the flux at 3.6 /im to the 
flux at 4.5 fiia at an age of 230 days correlates with the light curve decline rate at maximum brightness. Since 
the light curve shape at maximum is also correlated with the peak luminosity, we expect the CHI— CH2 color 
index at a fixed age is a good indicator of the radioactive nickel yield of a SNIa. 

• The series of optical spectra obtained with the LET show emission features fading at three different rates 
depending on the element and ionization state. Singly ionized iron-peak elements faded close to the ^^Co 
radioactive decay rate. Doubly ionized iron lines decayed at a rate similar the B-hernd fading. Finally, doubly 
ionized cobalt faded at more than twice the ^®Co radioactive decay rate. 

• The optical and near-IR spectra of SN 2011fe show large (> 1000 km s^^) velocity shifts in the singly ionized 
emission lines while doubly ionized lines showed significantly smaller offsets. 

It has long been known that the decline in radioactive energy input is the cause of the fading luminosity in the 
nebular phase. Here, we show that the doubly ionized species fade more quickly than the radioactive energy input 
and suggest that this additional fading is due to recombination. Additionally, doubly ionized cobalt sees a decline 



http: / /www. weizmann. ac.il/astrophysics/wisercp/ 
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in excitation, recombination losses and direct radioactive decay to iron and so it fades at the fastest rate of all the 
emission lines. The decay rate variation with element and ionization state demonstrated in the L BT spectra nicel y 
explains the range of observed decay rates seen in broadband filters during late-time studies of SNIa (|Lair "eFaIl[2006h . 

Our mid-IR Spitzer photometry is consistent with the slow fading of singly ionized iron in the CH2 band and doubly 
ionized iron or cobalt in the CHI filter. From the correlation between the IRAC color index with the radioactive nickel 
yield, we argue that Co III is probably the primary source of emission around 3.6 /im up to a year after explosion. The 
rapid fading of the Co III emission will eventually lead to Fe III or Fe II becoming dominant in CHI and a decelerating 
decay rate in the future. 

We al so measured the indiv idual line velocities of optical and NIR features and confirm the pattern of velocity offsets 
seen bv lMaeda et"all (|2010aD . That is, emissions from singly ionized species have higher velocity offsets than those 
from doubly ionized elements. According to Maeda et al. (2010b), this is because the singly ionized species reside in 
the dense inner region of the nebula where an offset deflagration occurred, while the doubly ionized atoms occupy the 
outer, low density and more symmetric regions of the nebula. Both the measures of the early and late spectra indicate 
SN 2011fe was a typical 'low velocity gradient' SNIa viewed at an angle of about 70° relative to the offset axis of the 
deflagration. 
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Table 1 

IRAC photometry for selected Type la Supernovae 



SN 


phase [days] 


CHI Flux [tijy] 


CH2 Flux 


2008Q 


222 


10.21 ± 1.08 


8.96 ± 1.51 


2008Q 


560 


0.38 ±0.28 




2008Q 


739 


< 0.42 




2009ig 


196 


13.25 ± 1.41 


2.44 ± 0.69 


2009ig 


375 


2.87 ± 1.05 


1.53 ±0.81 


2011fe 


146 


656.73 ± 13.60 


157.57 ±6.92 


2011fe 


167 


483.12 ± 11.91 


120.97 ±6.02 


2011fe 


234 


197.35 ± 7.99 


84.87 ±5.14 


2011fe 


351 


40.89 ±4.25 


34.28 ± 3.40 



Table 2 

Decline rates from least-squares fits to SN la IRAC photometry, and 
IR color differences. 



SN 


Ami5(B) 
(mag) 


A CHI 

[mag/lOO''] 


A CH2 

[mag/lOO"*] 


CH1-CH2 
@230 days [mag] 


2005df 
2008Q 
2009ig 
2011fe 


1.09 ±0.04 
1.36 ±0.06 
0.89 ±0.02 
1.07 ±0.06 


0.97 ± 0.03 
1.06 ±0.24 
0.93 ± 0.23 
1.48 ± 0.02 


0.65 ± 0.09 

0.28 ±0.36 
0.78 ±0.02 


-0.64 ±0.28 
0.34 ±0.22 
-1.14 ±0.41 
-0.44 ±0.08 



Table 3 

Nebular line identifications, velocities and decay rates. 







Identification'' 






ID Reference'^ 


Average Wavelength^ 


fncb 

(km s-i)= 


Decline Rate 
(mag/100<*) 




Fe HI] 4658 - 


f [Fe HI] 4701 4 


[Fe 


HI] 4734 


MIO 


4698 


-457±89 


1.41 




Fe HI] 4990 - 


f [Fe II] 5006 


B97 


4998 


-1638±180 


1.39 




Fe II] 5262 + 


[Fe III] 5270 






M10,B97 


5266 


-962±296 


1.55 




Co HI] 5890 


+ [Co HI] 5908 






B97 


5899 


-422±56 


2.01 




Co HI] 6129 


+ [Co HI] 6197 






B97 


6163 


+146±165 


2.33 




Co HI 6578 


+ [Fe II] 6456 + 


[Fe 


II] 6516 


B97,NIST 


6517 


+ 1391±1850'i 


1.65 




Fe II 


7155 + 


[Fe II] 7172 


M10,B97 


7154 


-1006±59 


0.85 




Ni II 


7378 + 


[Fe II] 7388 






M10,B97 


7383 


-1137±102 


0.81 




Fe II 


8617 






MIO 


8617 


-1559±101 


0.47 




Fe II 


1.257Aim 






MIO 


1.257/im 


-1223±95'= 






Fe II 


1.644/^m 






MIO 


1.644/im 


-1137±197<= 





^ Wavelengths arc given in A unless otherwise noted. 

^ References: (M10) [Ma"eda et al.|[20T0at (697) [Bowers et al.llT997l : (NIST'j iKramida et a'l]|2012l 

Veloc ities a re corrected for the MlOl redshift of z = 0.000804 listed by NED (NASA/IPAC Extragalactic Database; 
|http:/7ned .ipac.caltech.cdul 

The large uncertainty derives from an assumed 40 A error in fitting the highly-blended line profile. 
^ from LBT/LUCIFER spectrum obtained at age 234 days. 



